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Abstract 

We derive the high energy asymptotic behavior of gauge boson production cross 
section in a spontaneously broken SU (2) gauge theory in the next-to-next-to-leading 
logarithmic approximation. On the basis of this result we obtain the logarithmically 
enhanced two-loop electroweak corrections to the differential cross section of VP-pair 
production at ILC/CLIC up to the second power of the large logarithm. 

PACS numbers: 12.15.Lk 



1 Introduction 

The VP-pair production at e + e~ colliders plays a crucial role for testing the Standard Model 
of electroweak interactions. At LEP2 this process has been used for the determination of 
the VP-boson mass Mw, a fundamental parameter of the standard model, through PP-boson 
reconstruction with an uncertainty of 40 MeV pQ. Furthermore, the triple gauge boson cou- 
pling as predicted by the non-Abelian gauge theory has been verified within a few percent. 
The experimental study of the VP-pair production at the International Linear Collider (ILC) 
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is expected to improve the accuracy of the mass determination to 7 MeV due to much higher 
luminosity [2j. Moreover, the advent of ILC will give access to the new high energy domain 
where the cross section is increasingly sensitive to the triple gauge boson coupling and W- 
pair production could be used as a probe of the non-Abelian structure of the electroweak 
interactions and of possible gauge boson anomalous couplings. To match the experimental 
accuracy, the theoretical analysis has to take into account the electroweak radiative correc- 
tions. The one-loop corrections to the cross section of the on-shell IV-pair production have 
been evaluated by different groups [3|B|5|l6] already decades ago. The calculation of the one- 
loop corrections to the W^-boson mediated e + e~ — > 4/ processes has been performed in the 
double pole approximation in Ref. [7] and incorporated into the event generators YFSWW 
[8] and Racoon WW [9j. Recently the full analysis has been completed [TUj. These results 
ensure an accuracy significantly better than one percent when the characteristic energy \J~s 
of the process is about the gauge boson mass. However, once y/s is far larger than M w , 
the cross section receives virtual corrections enhanced by powers of electroweak "Sudakov" 
logarithm hx{s/M^ z ), which at the energies of about one TeV have to be controlled to 
two loops in order to keep the theoretical error below one percent. This is even more valid 
for energies of 3 TeV anticipated for the CLIC project [llj. In the case of light fermion 
pair production these corrections are already available through the next-to-next-to-next-to- 
leading logarithmic (N 3 LL) approximation, i. e. including all the two-loop logarithmic terms 
[12)13)14)15)16] . For e + e~ — > W + W~ production only the leading-logarithmic (LL) and the 
next-to-leading logarithms (NLL) are known so far [l~2")17)18Jfl~9] . 

In the present paper we extend the analysis of H^-pair production to the next-to-next-to- 
leading logarithmic (NNLL) approximation following the approach developed in Refs. [Lf)14M15|[T6i 
for the four-fermion processes. The limit of the small-angle production, which could be in- 
teresting in the case of the transverse gauge bosons because the corresponding cross section 
is peaked in the forward direction, remains beyond the scope of the present paper. In the 
next section we outline the approach and derive the NNLL corrections to the differential 
cross section of the gauge boson pair production in a spontaneously broken SU(2) model 
which emulates the massive gauge boson sector of the Standard Model of electroweak inter- 
actions. We generalize the result to the £77(2) x £7(1) gauge theory with a heavy top quark 
in Section [3J A brief summary and conclusions are given in Section HI 

2 High energy asymptotic of the massive gauge boson 
production 

Let us consider as a toy model the spontaneously broken SU(2) gauge theory with the Higgs 
mechanism for gauge boson mass generation and with massless left-handed fermion doublets. 
The model retains the main features of the massive gauge boson sector of the Standard Model. 
In this case the result can be presented in a simple analytical form and constitutes the basis 
for the further extension to the full electroweak theory. We study the process of gauge boson 
pair production in fermion-antifermion annihilation at high energy and fixed angle when all 



2 



► [WWW 

— * — k/w/w 

(a) 

Figure 1: The diagrams represent (a) transverse and (b) longitudinal gauge boson pair production 
in fermion-antifermion annihilation at high energy in the Born approximation. 

the kinematical invariants are of the same order and far larger than the gauge boson mass, 
\s\ ~ \t\ ~ \u\ ^> M 2 . In this limit the asymptotic energy dependence of the field amplitudes 
is dominated by Sudakov logarithms [2"Uf2"T] and governed by the evolution equations. The 
method of the evolution equations in the context of the electroweak corrections is described in 
detail for the fermion pair production in Ref. [16]. The derivation of the evolution equations 
[22|23f24j applies to any process of wide-angle production or scattering of on-shell particles 
when the characteristic momentum scale is far larger than the mass scale. It is entirely 
based on (i) the properties of hard momentum region and (ii) ultraviolet renormalization 
of the light-cone Wilson loops. Therefore it depends neither on specific infrared structure 
of the model nor on the specific choice of the external particles (for the extension to the 
processes with arbitrary number of external particles see Ref. [25]). Thus, the approach of 
Ref. [16] directly extends to the gauge boson production as briefly described below. The 
only potential subtle point in the analysis of gauge boson production is that the effects of 
spontaneous symmetry breaking can change the asymptotic states as it happens with photon 
and Z-boson in the standard model. This would require an additional consideration. We, 
however, restrict the analysis to the production of V^-bosons which have the same gauge 
quantum numbers in broken and symmetric phases and do not encounter this problem. 

Due to helicity conservation a pair of either transverse or longitudinal gauge bosons can be 
produced in the high energy limit. The transverse gauge bosons behave like vector particles 
in the adjoint representation while the longitudinal gauge bosons, as a consequence of the 
equivalence theorem, behave like scalar particles in the fundamental representation. The 
structure of the Sudakov logarithms in these cases is significantly different and we consider 
them separately. 

2.1 Transverse polarization 

The transverse gauge bosons are the true vector particles and the Born amplitude in this case 
is given by the t-channel and w-channel fermion exchange diagrams, Fig. [Tfr. It is convenient 
to introduce the functions Z^ and Za which describe the asymptotic dependence on the 
large momentum transfer Q of the fermion scattering amplitude in an external singlet vector 
field and of the vector boson in an external scalar singlet field, i.e. of the respective form 
factors in the Euclidean region. In leading order in M 2 /Q 2 these functions are known to 
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satisfy the following linear evolution equation [22|23ll2^ 



a 



d\nQ 2 
with the solution 



Z; 



Q 2 dx 



M 2 X 



iMx)) + 0(«(Q 2 )) + Ua{M 2 )) z % 



Q 2 dx 




dx' 

— iMx')) + (i(a(x)) + 6(«(M 2 )) 

M 2 T 



(2) 



which satisfies the initial condition Zi\q 2=M 2 = 1- Here the perturbative functions 7»(at) etc. 

are given by the series in the coupling constant a(jj 2 ), e.g. ji(a) = J2%Li( a l '^) n Ji ■ For 
the amplitude of transverse boson production At let us introduce the reduced amplitude At 
so that 

At = Z^Z a At . (3) 

Due to the factorization property of the Sudakov logarithms associated with the collinear 
divergences of the massless theory [26J the reduced amplitude satisfies the simple renormal- 
ization group like equation [27ll28f29] 
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;A T = Xt{<Q 2 ))At 
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where Q 2 = —s and \ T is the soft anomalous dimension matrix acting in the space of the 
isospin amplitudes. The solution of the above equation is given by the path-ordered exponent 



At = Pexp 



Q 2 dx 



M 2 X 



Xr(a(z)) 



At(«(M 2 )) 
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where Atq determines the initial conditions for the evolution equation at Q = M. By 
calculating the functions entering the evolution equations order by order in a one gets the 
logarithmic approximations for the amplitude. For example, the LL approximation includes 
all the terms of the form a n C 2n and is determined by the one- loop coefficients 7 l - 1 ' ) . The NLL 
approximation includes all the terms of the form a n £ 2n ~ m with m = 0, 1 and requires the 
one-loop coefficients Q , Q , and as we U as the one-loop running of a in ji(a), and so 
on. To get the NNLL terms a n C 2n ~ 2 one needs in addition the two- loop coefficient 'y^, the 
two- loop running of a in 7»(a) and the one- loop contribution to Ato- 

The anomalous dimensions 7(a), ((a) and x( a ) are mass-independent and can be associ- 
ated with the infrared divergences of the massless (unbroken) theory. From the QCD result 
(see e.g. [20] and references therein) adopted for the specific case of SU(2) gauge group, rif 
chiral quarks, and one scalar in the fundamental representation we get 
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and 7^ = 87I /3. The matrix x^P can be extracted from the results of Refs. [3T1I32 ]. In 
the isospin basis (er a cr 6 , cr b (T a , 5 ab ■ 1) it takes the form 



X? 



/-2(ln(x_)+m) ln(f±)) 

-2(ln(x + )+m) ln(f^ 

V (ln(x+) + in) (In(x_)+Z7r) 



(7) 



where x± = (1 ± cos#)/2 and 9 is the production angle. Note that in this basis the Born 
amplitude up to a common factor is given by the vector l/x+,0). At the same time 

the functions £i(a) and At(o) do depend on the infrared structure of the model and require 

and from the 
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the calculation in the spontaneously broken phase. For example £, 

result of Ref. [6] we obtain $ = 0- To emulate the e + e~ — > W^W^ process one has to 
project the amplitude on the relevant initial and final isospin states, which is straightforward. 
The Born cross section in the high energy limit reads 



~dn 



a 2 (M 2 ) x+(x\ + x 2 _) 



As 



and is peaked in the forward direction. We define the perturbative expansion for the differ- 
ential cross section in the MS renormalized coupling constant a = a (M 2 ) as follows 
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Expanding the Sudakov exponents to NNLL order for the one- and two-loop corrections we 
get 
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where (3q = 43/6 — H.//3 is the one-loop beta-function, Qj = Refxy' 1 ]^, C(s) = ln(s/M 2 ), 
and 5qt is the nonlogarithmic part of the one-loop corrections which can be extracted from 
the result of Ref . [6] . For the Higgs boson mass Mjj — M it takes a simple form 
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Substituting the values of the coefficients for n/ = 12 we obtain 
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Note that in contrast to the four-fermion processes, the cross section of the gauge boson 
production depends on the Higgs boson mass already in the NNLL approximation. 



2.2 Longitudinal polarization 

The equivalence theorem relates the amplitude of the longitudinal gauge boson production 
e + e~ — > W£W£ to the production of the Goldstone bosons e + e~ — > 4> + (fi~- The Born 
amplitude is now given by the s-channel annihilation diagram, Fig. [TJ. The analysis of 
the high energy asymptotic for the longitudinal gauge boson production goes along the line 
described in the previous section and is very similar to the one for fermion pair production 
|14j . Instead of Z& one should use the function Z^ which correspond to the scalar particle 
scattering in an external singlet vector field. The necessary parameters of the evolution 
equation read 7^ = 7^ , Q/ = 3, and from the result of Ref. |6J we get Q 1 -* = 0. The 
structure of the reduced amplitude is also different. In the isospin basis (<x a ® cr a , 1 ® 1) the 
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one-loop matrix of the soft anomalous dimensions takes the form 

/-4(k( I+ )+ix) + 21„(|j) 41n(|±)\ 

XL ~[ iHt) o )' ( ) 

Note that in this basis the Born amplitude is given by the vector (1, 0). The corresponding 
Born cross section reads 

i R 9 / \ 



daf a 2 (s) x + x. 



(\A> l.s I 

It has a maximum at 9 = 90°. We proceed as in the case of transverse polarization and 
obtain the one- and two-loop NNLL corrections to the differential cross section 

Sf = {if + if) A*) + 2 [Cf + Cf + if + & + 1$ + <] C{s) + *f (17) 

*f = \ {if + lf) 2 C\s) + 2 [(f + Cf + if + if + & + 4« } - \Po] {if + if) 

*£\s) + [if + if + 2 {cf + cf + if + if) {cf + cf + 4 1} + & 

+2/ff + 81$) - P + < + Cf + Cf) + If' + 1$ (4/!i } + 1$ + 4/g) 

(/Si ) + 44i ) ) 2 + 4 ) (7i 1) + 7i 1) )]^(^) 
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where if = Re[xf)ij- From the result of Ref. [6] for Mh = M we obtain^ 
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For nf = 12 this gives 
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ln 2 (x + ) — 8 ln(x_) ln(x_| 
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(21) 



1 The equivalence theorem holds up to the field renormalization factor (see e.g. Ref. [33]) which affects 
the initial conditions for the evolution equation. Thus one has to use the explicit result for the longitudinal 
W-boson production rather than the equivalence theorem to get the momentum and angular independent 
term in . 
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3 VF-pair production in e + e annihilation 

In the standard electroweak model the perturbative expansion involves the SUl{2) coupling 
constant a ew and the U(l) hypercharge coupling constant ay. We eliminate the latter by 
means of the relation ay = tan 2 6y/ a ew , where 9\y is the weak mixing angle, and consider 
the one-parameter series for the cross section in a ew of the form of Eq. ([9l). In the high 
energy limit the transverse gauge bosons are produced only in annihilation of the left-handed 
electron-positron pair. The corresponding Born cross section is given by Eq. (JHJ) with a 
replaced by a ew . The longitudinal gauge bosons can be produced in the annihilation of the 
electron-positron pair of both chiralities. In the case of the left-handed initial state fermions 
the Born cross section gets the contribution from the SUl{2) and the hypercharge virtual 
gauge bosons and reads 

d(7 -L = 1 al w (s)x+x- , . 

dtt cos 4 #vk As A ' 1 ' 

For the right-handed initial state fermions the Born cross section is saturated by the hyper- 
charge gauge boson 

d °?L _ „ .4 - clff B 



dn w »"-3jr (23) 

The analysis of the radiative corrections in the Standard Model is complicated by the presence 
of the mass gap and mixing in the gauge sector and by the large top quark Yukawa coupling. 
In next section we use the method of Ref. [T2~yT6] to separate the electroweak and QED 
Sudakov logarithms. In Sect. 13.21 we extend the evolution equation approach to the Yukawa 
enhanced contributions. In Sect. 13.31 we present the final numerical results for the two-loop 
NNLL corrections to the differential cross sections. 



3.1 Separating QED Sudakov logarithms 

The electroweak Standard Model with the spontaneously broken SUl(2) x U(l) gauge group 
involves both the massive W and Z bosons and the massless photon. The corrections to 
the fully exclusive cross sections due to the virtual photon exchange are infrared divergent 
and should be combined with soft real photon emission to obtain infrared finite physical 
observables. We regularize the infrared divergences by giving the photon a small mass A. 
Thus besides the electroweak Sudakov logarithms discussed above the radiative corrections 
contain the QED Sudakov logarithms of the form ln(Q 2 /A 2 ). To disentangle the electroweak 
and QED logarithms we use the approach of Ref. [L2fl4lll6| . While the dependence of the 
amplitudes on the large momentum transfer is governed by the hard evolution equations 
(c.f. Eqs. ([Hill)), their dependence on the photon mass is governed by the infrared evolution 
equations [12]. In the limit A 2 <C Myy, in 2 ^ Q 2 the infrared evolution equations in the 
full theory are the same as in QED and the solution to the NNLL accuracy in the massless 
fermion approximation mj — (/ 7^ t) is given by the factor 
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27 7T 

+0(al)j, (24) 

where a e is the MS QED coupling constant. The NNLL approximation for U can be ob- 
tained from the result for the fermion-antifermion production [TJ] by proper modification 
of the QED anomalous dimensions. It is convenient to normalize the QED factor so that 
=A 2 =M 2 =1. In order to cancel the singular dependence on the photon mass, the 
QED Sudakov exponent ( 1241) should be combined with the the real photon emission, which 
is also of pure QED nature if the energy of emitted photons is much smaller than My/- 

Two sets of equations completely determine the dependence of the amplitudes on two 
dimensionless variables Q/Mw and Qj A up to the initial conditions which are fixed through 
the matching procedure. To get the purely weak logarithms one subtracts the QED expo- 
nent fT2~4|) from the exponent given by the solution of the hard evolution equation. This can 
naturally be formulated in terms of the functions parameterizing the solution. The functions 
7, (, and x are mass- independent. Therefore the anomalous dimensions parameterizing the 
purely weak logarithms can be obtained by subtracting the QED contribution from the re- 
sult of the unbroken symmetry phase calculation to all orders in the coupling constants. In 
the order of interest they can be found in or easily derived from the result of Ref . [T4] . For 
example, we have 



(2) (2) I 8UU • 2n /q rr \ 

Ta = Ta \su(2) ~ -^f 8m w, (25) 
if = lf ) \su(2) + ytan 2 ^-^sin 2 ^, (26) 



and so on. Here the SU(2) contributions are given by the results of Sect. [2] with nj = 12. 
The only new ingredient in comparison with the light fermion pair production [TJ] is the 
effect of the large Yukawa coupling of the third generation quarks on the longitudinal gauge 
boson production which is considered in the next section. 

On the other hand the functions £ and Ao are infrared sensitive and require the use of 
the true mass eigenstates of the Standard Model in the perturbative calculation. In the 
NNLL approximation one needs the one-loop contribution to these quantities which can be 
extracted from the result of Ref. [6]. For example, for the left-handed initial state fermions 
we find 



(1) (1) I -4cOS 2 g l y + 8cOS 4 ^ / M 2 



w , 



4" +4" - (i - 2 t; m M #i. 

COS" 1 W \ Myy 



At the same time the Aq term results in the one-loop nonlogarithmic contribution to the 
cross section. The corresponding expression directly follows from Ref. [6] and is rather 
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Figure 2: The one- loop diagrams contributing to the anomalous dimension matrix £. The arrow 
lines correspond to the third generation quarks. The dashed lines correspond to the Higgs, neutral 
or charged Goldstone bosons. The black square represent an external singlet vector field 



lengthy so we do not give it explicitly. Note that in Ref. [6] the result is presented in the 
on-shell renormalization scheme and in the limit A <C m e . We convert it to MS scheme and 
to the massless electron approximation using the formulae of Refs. [Mf35] . 

3.2 Top quark Yukawa coupling effects 

The large Yukawa coupling of the third generation quarks to the scalar (Higgs and Goldstone) 
bosons results in specific logarithmic corrections proportional to m^/M^. The high energy 
evolution of the form factors in a theory with Yukawa interaction is completely analogous 
to the one of 3 scalar theory in six dimensions, see the second paper of Ref. [23]. The 
structure of factorization and evolution equations is much simpler than in a gauge theory 
because Yukawa interaction itself does not contribute to the anomalous dimension 7^(0;) and 
results only in single logarithmic corrections completely determined by the ultraviolet field 
renormalization of the external on-shell particles. These corrections can be taken into account 
through the modification of the evolution equations for the corresponding ^-functions. The 
analysis is straightforward but complicated because the Yukawa interaction mixes evolution 
of the quark and scalar boson form factors and in general does not commute with the 
SU(2) and hypercharge couplings. However, due to the factorization of the double Sudakov 
logarithms, the Yukawa enhanced contribution to NLL approximation is given simply by 
the product of the one-loop Yukawa corrections and the double logarithmic exponent as 
observed in Ref. [36]. The structure of the NNLL contribution is much more complicated 
and we restrict the analysis to a simplified model with sin8w = 0, i.e. with no hypercharge 
interaction. Let us introduce the following five-component vector in the space of ^-functions 
Z = (Z^, Z x , Z b _, Z t _, Z t+ ), where the subscript + (— ) stand for the right (left) quark 
fields and Z x corresponds to the transition of the Higgs boson into the neutral Goldstone 
boson in the external singlet vector field. The evolution equation for this vector takes the 
form 



«91nQ 2 




(x)) + C(a ew (Q 2 ), atYuk{Q 2 )) + £(a e «,(M^)) Z , 



(28) 



with the solution 



„ ^ f f Q2 dx \ f x dx 
Z = Pexp I — / — 
[Jm^ x [Jm^ x' 



T j(a ew (x')) + C(a ew (x), a Yu k(x)) + ^(a ew (M^)) \ Z , 



(29) 
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where = (—3/2) • 1 and £ = 0. The anomalous dimension matrix £ includes all the 
dependence on the Yukawa coupling ay u k- We eliminate the latter by means of the relation 



&Yuk 



2M, 



h-otew, and consider the one-parameter series for the anomalous dimension in a € 



The one-loop coefficient reads 



/12 0\ 

12 
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9 

V 0/ 



+ 



mi 



4M 2 , 



/ 6 -6\ 

6 -6 

1 -1 
10 0-1 

V-l -1 -1-1 / 



(30) 



where the first term representing the pure SUl(2) contribution follows from the result of 
Sect. 12.11 and the second term represents the Yukawa contribution. It can be extracted 
from the known one-loop result (see e.g. Ref. |36,37j). The relevant diagrams are given 
in Fig. [2j Note that instead of the ^-functions associated with the form factors one can 
directly consider the ultraviolet field renormalization. In this case the non-diagonal form of 
the anomalous dimension matrix is due to the mixing of the bilinear quark and scalar boson 
operators, which is specific for Yukawa interaction and is absent in a gauge theory. 

The first two diagrams correspond to the mixing of the quark and the scalar boson form 
factors. Moreover the Yukawa coupling changes quark chirality and/or flavor and the last 
diagram in Fig. [2] corresponds to the pure mixing of Z^, Z t _ an d Z t+ functions. As a 
consequence, all the diagonal matrix elements in the second term of Eq. (I3"U|) vanish. 

The proper initial condition for the evolution equation which corresponds to the Born 
amplitudes of the quark and scalar boson production in e + e~ annihilation is given by the 
vector Zq = (1, —1, —1, 1,0). In NNLL approximation one needs also the one-loop running 
of the Yukawa coupling with the corresponding beta- function /3^ uk = | — j^h- By expanding 
the solution for the component Z^ we obtain the two-loop corrections enhanced by the second 
or fourth power of the top quark mass. Note that in the production amplitude one has to 
take into account also the interference between the one-loop Yukawa contribution to Z^ and 
the one-loop logarithmic term in the reduced amplitude and the electron Z^ function. The 
two-loop NNLL Yukawa enhanced contribution to the cross section reads 
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(2) 

NNLL 



Yuk 



3 m\ 
2M^ 



+ — f (301n(x_)-61n(x+)-27) 



Ml 



C 2 (s). 



(31) 



This expression approximates the full result up to the terms suppressed by sin 2 9w ~ 0.2. 



3.3 Numerical results 

We adopt the following input values j3S] M = M w = 80.41 GeV, M H = 117 GeV, m t = 
172.7 GeV for the masses and sin 2 9w = 0.231, a ew = 3.38- 10~ 2 for the MS coupling constants 
renormalized at the scale of the gauge boson mass. Note that the coupling constants in the 
Born cross section of the longitudinal gauge boson production are renormalized at the scale 
>/5. 
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Figure 3: The one- loop logarithmic corrections to the differential cross section relative to the Born 
approximation at yfs = 1 TeV as functions of the production angle for (a) transverse and (b) 
longitudinal polarization of the gauge bosons. 



Transverse polarization. We obtain the following one and two-loop NNLL corrections 
to the cross section 



4 1} 



4 2) 



-4.73 C 2 (s) 4 

4.70 2.35 + 1.95z 



6.15 + 4.00 — j ln(x_) + 2.15 ln(a; + ) + 4.43 



C(s) 



+ 



+ 



X _i 



xt + x\ 



\yy\x-) + 



4.00 



X j 



ln(x_) ln(x + ) + 



3.00 x_ , ,, 

ry*^ I /yi2 



In (xa 



4.95 - 9.65 x_ \ , . . nrA1 . . 2.35x+ 
+ | 0.54 + — — 5 I ln(s_) - 0.54 ln(x+) - — \ - 0.82 , 



,->-> 2 i ,->t2 
x r x_|_ 



x_ + an 



(32) 



1891x _\ 

29.08 : ln(ar_) - 10.17 ln(a; + ) - 14.62 

22.21- 12.61 rr_ 4.00x 2 11.11 + 9.22x_\ , 
18.92 + + 5 5 I ln 2 (a;_ 



11.17£ 4 (s) + 
+ 

- ( 9.24 + 



£ 3 ( S ) 



18.91 3.39x_ 



X2 I 2 

— \ Xl 



/ _ 11.40 x_ 23.40 - 45.61 xS 
- 15.26 h - 



, \ i / \ / 14.18 x_\ , 2/ 
ln(x_) m(> + ) + 2.32 * - ) \n 2 (x A 



x 2 _ + xl 



Xa 



x 2 _ + x\ 



ln(x_) + 3.87 ln(o: + )+ 1 ^ 11X 2 + 



-3.60 



C 2 (s) 



(33) 



The one-loop result is well known j6] and is given here to show the structure of the logarithmic 
expansion, see Figs. [3fr and Hh. In Figs. \&p> and Eh the values of different logarithmic two- 
loop contributions as well as their sum are plotted as functions of the production angle 
at the center of mass energy of 1 TeV and 3 TeV, respectively. The two-loop subleading 
contributions exceed the LL one in absolute value in the small angle region. However, due to 
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the partial cancellation between the NLL and NNLL terms the total NNLL approximation 
is close to the LL one. It has a fairly flat angular dependence and amount to about 5% for 
y/s = 1 TeV and 15% for y/S = 3 TeV. 

Longitudinal polarization. For the left-handed initial state fermions the one and two- 
loop NNLL corrections to the cross section read 



f( 2 ) 



-2.38 C 2 (s) + [-6.91 ln(ar_) + 0.75 ln(ac+) 

9 1 Q fl f>5 

-±±1 ]^( x _) + l n 2( x+ ) + 0.19 (ln(x_ 



-•4 



3.48] C(s) 
-ln(x+)) + 36.85. 



(34) 



2.82 £ 4 (s) + [16.41 ln(x_) - 1.79 ln(x + ) + 11.87] £ J (s) 



+ 



5.20 \ , „ 
18.38 + ln 2 {x. 

x + J 



0.28 - — I ln 2 (x. 



x. 



+3.11 ln(x_) ln(a; + ) + 49.(10.) ln(ar_) - 18.(4.) \n(x+) - 128.(20.) 



C 2 (s) . (35) 



In the two-loop NNLL contribution the error bars indicate the uncertainty due to our approx- 
imation of the Yukawa enhanced contribution. The structure of the logarithmic corrections 
differs from the case of the transverse polarization as one can see on Figs. [3fr, Hfc and0>. The 
sum of the two-loop logarithmic terms is strongly angular dependent and varies between -3% 
and 2% for <fe = 1 TeV and between -7% and 8% for = 3 TeV. 

For the right-handed initial state fermions the Born cross section is suppressed by the 
factor 4sin 4 #yi/ « 0.2 in comparison to the left-handed case. Moreover the two-loop loga- 
rithmic corrections turned out to be about 3 • 10~ 3 for all the scattering angles. Thus they 
are of no phenomenological importance and are not presented here. 
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Figure 5: The two- loop logarithmic corrections to the differential cross section relative to the 
Born approximation at yfs = 1 TeV as functions of the production angle for (a) transverse and ib) 
longitudinal polarization of the gauge bosons. 



4 Summary 

In the present paper we employed the evolution equation approach to analyze the high energy 
asymptotic behavior of the gauge boson production through the annihilation of fermion- 
antifermion pair in the spontaneously broken SU(2) gauge model. The result has been used 
to compute the two-loop NNLL electroweak corrections to the differential cross section of Im- 
pair production in e + e~ annihilation. The corrections are comparable and even exceed the LL 
terms depending on the production angle. The structure of the corrections is different for the 
transverse and longitudinal boson production. In the first case we observe the cancellation 
between the huge NLL and NNLL contributions so that the sum is dominated by the LL 
term and amounts of about 5% at ^/s ~ 1 TeV and 15% at */s ~ 3 TeV. For the longitudinal 
bosons the corrections exhibit significant cancellation between the LL, NLL and NNLL terms 
so that the sum does not exceed 2% in absolute value for <Js ~ 1 TeV. The cancellation 
becomes less pronounced at higher energy. The uncertainty of the theoretical prediction for 
the on-shell 14^-pair production at ILC is now determined by the unknown two-loop linear 
logarithmic terms. For the fermion pair production such terms are know to contribute about 
1-2% of the cross section [16]. This value can be used as a rough estimate of the accuracy 
of our approximation. We should emphasize that our approximation breaks down at small 
production angles where the Regge logarithms ln(— t/s) become large. 
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